Helical edge states coupled to a spin bath: Current-induced magnetization 
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We study current carrying helical edge states in a two-dimensional topological insulator coupled 
to an environment of localized spins, i.e. a spin bath. The localized spins mediate elastic spin-flip 
scattering between the helical edge states, and we show how this induces a spin-bath magnetization 
for a finite current through the edge states. The magnetization appears near the boundaries of 
the topological insulator, while the bulk remains unmagnetized, and it reaches its maximal value 
in the high bias regime. Furthermore, the helical edge states remain ballistic in steady state, if no 
additional spin-flip mechanisms for the localized spins are present. However, we demonstrate that 
if such mechanisms are allowed, then these will induce a finite current decrease from the ballistic 
value. 

PACS numbers: 72.10.Fk, 72.25.-b, 73.63.-b, 75.76.+j 



I. INTRODUCTION 

Gapless helical edge states exist at the boundary of 
a two-dimensional (2D) topological insulator (TI)i~— . 
A pair of helical edge states are counter propagating, 
and different edge states with opposite wave numbers, 
k and —k, constitute a Kramers pair. The prime ex- 
ample of a 2D TI has been recently realized experimen- 
tally by the group of Molenkamp-t~- in HgTe quantum 
wells. Evidence of edge states in both two-terminal^ and 
non-local^ transport measurements were found and, fur- 
thermore, the relation between the spin of the helical 
edge state (HES) and propagation direction was exper- 
imentally established^.. The existence of the TI state 
in a HgTe quantum well beyond a certain critical well 
thickness was predicted by Bernevig, Hughes and Zhang 
(BHZ) 8 by constructing a minimal model - similar to a 
massive Dirac model - describing the basic physics. 

The fact that the HESs come in Kramers pairs means 
that elastic scattering from wave-vector k in one HES to 
— k in the other HES within a pair cannot be induced by 
time-reversal invariant potentials. Therefore scattering 
by e.g. impurities between a pair of HESs is strongly sup- 
pressed, which results in ballistic transport and quantized 
conductance of e 2 /h per pair of HESs as observed^. In- 
elastic scattering mechanisms*^— and scattering mech- 
anisms breaking time-reversal invariance can, however, 
induce scattering between the HESs. For instance, the 
effect of a single magnetic impurity on the conductance 
through a pair of HESs have been considere d 12 ^ 3 . In 
this case, Tanaka et alM- showed that even though a sin- 
gle magnetic impurity breaks time-reversal invariance, it 
does not have any effect on the dc conductance. Also the 
RKKY interaction mediated by HESs has been studied^. 



A. Qualitative considerations on the current 
induced magnetization and the current change 

In this paper, we consider a 2D TI coupled to an envi- 
ronment of localized spins, i.e. a spin bath. We focus on 




FIG. 1. (color online) (a) A pair of helical edge states coupled 
to a spin bath (blue arrows) between two leads with chemical 
potentials /i£ and /j,r, respectively. The pair of helical edge 
states exists at the interface between a 2D topological insula- 
tor (white region) and an ordinary insulator (grey area), (b) 
and (c) show the dominant spin-flip scattering processes be- 
tween the edge states for fit > hr and (j,r > /il, respectively, 
for an unmagnetized bath. Here the right- (left-) moving edge 
state is occupied up to fiL (fJ-n), as indicated by the thick blue 
(red) lines. Since the inter edge state scattering is mediated 
by a flip of a localized spin in the bath, it will lead to an 
excess number of spin up in the bath for [il > (ir in (b) 
(and vice versa for /j,r > hl in (c)), and in terms give a finite 
magnetization. 



a single pair of HESs as seen in Fig. [TJa), even though 
a real device has a pair of HESs at each boundary. This 
can be done without loss of generality as long as the 
boundaries are not close together (as for instance in a 
point contact geometry^). Furthermore, for simplicity 
we consider the case of a spin-1/2 bath, which is not 
essential for the physics discussed in this paper. 

The spin bath breaks time-reversal symmetry (from 
the point of view of the carriers in the HESs) and there- 
fore enables the possibility of elastic scattering between 
the HESs. We argue that current through the HESs will 
induce a magnetization in the spin bath near the bound- 
ary region of the 2D TI. The magnetization comes about 
due to angular momentum conserving scattering between 
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the HESs. The occupation of the HES is determined by 
the contact, where it originates. This means that the 
right moving HES is occupied up to the chemical poten- 
tial hl of the left contact and vice versa (see Fig. [I}. 
Therefore, a finite bias voltage, say (1l > A'-Ri opens an 
energy window favoring scattering from the right-moving 
spin up HES to the left-moving spin down HES, see 
Fig. [TJb). This scattering is mediated by a spin flip in 
the bath in the opposite direction, I to j, and therefore 
this rate of inter-HES scattering depends on the num- 
ber of spin down in the bath. The scattering between 
the HES will therefore dynamically change the number 
of spin down and up in the spin bath until a steady state 
is reached. In other words, the magnetization in the spin 
bath builds up to compensate the increased phase space 
for the inter-HES scattering due to the bias. In particu- 
lar, the spin bath can magnetize completely in the high 
bias regime, where the phase space for one of the two 
inter-HES scattering processes is suppressed completely. 
Since the HESs only exist on the boundary, the spins lo- 
calized in the bulk of the 2D TI are not affected by the 
current through the HESs. 

The spin-flip scattering between the HESs also change 
the propagation direction of the carrier, i.e. it is a 
backscattering process. However, if no additional spin- 
flip mechanisms are present for the localized spins, then 
the HESs remain ballistic in steady state, because once 
a spin in the bath has mediated a transition between the 
edge states by flipping from -say- down to up, then it 
cannot mediate another transition. Nevertheless, if ad- 
ditional spin-flip mechanisms are feasible to randomize 
the direction of the localized spins, then this will induce 
a finite steady state current change, since backscattering 
between the edge states now will try to compensate the 
randomizing of the spins. 

One possible experimental realization of the spin bath 
is magnetic impurities, e.g. Mn ions in a HgTc quantum 
welli^. In this case, low concentration of magnetic im- 
purities is required not to hinder the existence of HESs. 
However, the physics and phenomenon discussed here is 
of a rather generic nature for any 2D TI coupled to an 
environment of localized spins. 

II. THE SPIN-BATH DYNAMICS AND THE 
CURRENT 

Next, we describe in detail how the spin bath and the 
transport through the HESs are connected. Without the 
spin bath, the electric current through a single pair of 
HESs is ballistic such that 

/ (0) = ^(ML-M«), (1) 

where (/ir) is the chemical potential of the left (right) 
lead and e > is the elementary charge. 

The spin-bath mediated scattering between the HESs 
is a backscattering of a single particle, since the two HESs 



are countcrpropagating. The electric current change 51 
due to the spin bath is given by the rate of change in the 
number of left (or right) movers. This is, in terms, given 
by the rates T a ^ a i for transferring a particle from the 
HES a 1 to the HES a, i.e. 

S/=(-e)(IV^-IV t ), (2) 

and the total current is I = I 1 - ) +SI. The scattering rates 
Tu^o-i depend on the magnetization of the spin bath. For 
instance, the more localized spins with spin up, the larger 
the rate r-|-<_^ and vice versa. (Detailed expressions are 
given below, see e.g. Eq.®.) 

Now we turn to the dynamics of the spin bath. Every 
time one particle is scattered between the HESs, a single 
localized spin is flipped. This means that the number 
of each spin species in the bath covered by the HES, 
N a , will change dynamically with the same rates T a <- a i 
as the inter-HES scattering. Therefore, the time evolu- 
tion of the number of localized spin species a, N a , due 
to the scattering between the HESs can be described by 
the rate equations dtN^ = — r^j., and dtN± = 

r^<_^ — r^-)-, where the total number of localized spins 
covered by the HESs, N s = N t + N±, is fixed. (Note that 
the rate indices refer to the HES spin flip, which is oppo- 
site to the spin flip in the bath.) For convenience, we nor- 
malize the magnetization such that the maximal (mini- 
mal) magnetization, where all the localized spins are in 
the up (down) state, is 1 ( — 1). Therefore, we write the 
magnetization as M. = (N-f — N±)/N s , which is also often 
called magnetic polarization. Therefore, subtracting the 
two rate equations for iV-j- and N±, the rate equation for 
the magnetization becomes dtM. = 2(rj_^ — T-j-^) /N s . 
However, this is considering only the spin flip in the 
bath stemming from the scattering between the HESs. 
Other - presumably much weaker - mechanisms might 
also flip the localized spins such as dipole-dipole inter- 
actions within the spin bath and spin-phonon coupling 
e.g. relevant for Mn ion a 17 ' 18 . Such mechanisms will try 
to equilibrate the number of spin up and down in the 
bath, and thus drive the magnetization towards zero. We 
include this in the time evolution of the magnetization 
by a phenomenological term —T r M /N s similar to the 
relaxation-time approximation^, i.e. 

dtM = J-(iVt - iVi) - j^r r M. (3) 

The phenomenological term is divided by N s such that 
T r is a spin-flip rate per localized spin and thereby com- 
parable to T . 

Therefore, it is now evident that for T r = 0, these 
simple rate equations lead to 5I(t) — ^N s dtM, such that 
in steady state, dtM = 0, there is no current change, 

61 = (for T r = 0), (4) 

and the HESs remain ballistic. Physically, the magne- 
tization builds up to compensate the difference in scat- 
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tcring rates between the two HESs (i.e. — is 

required in steady state for T r =0). 

Taking the additional weak spin-flip mechanisms in the 
bath into account, T r ^ 0, the current change is no longer 
zero in the steady state, but found to be 



SI = -j^rM + 



(5) 



by inserting the current change (0) into dtM = 0. Phys- 
ically, the current change is a result of the competition 
between the additional spin flip mechanisms within the 
bath and the spin flips due to the inter-HES scattering. 
Below, the scattering rates are found such that the mag- 
netization and current change can be studied in greater 
detail. 



III. THE HELICAL EDGE STATES AND THEIR 
COUPLING TO A SPIN BATH 



We model the HESs by the eigcnstates 
( p kt (x,y) = -^=e ik *f k (y)\V, 



(6a) 
(6b) 



where f k (y) is the (real) transverse wavefunction of width 
W y localized at the boundary of the 2D TI and L is the 
length of the HES, sec FigfTJa). The energies are e ka = 
Ea+shvok, where s = +1(— 1) is for spin a =t (|), v n > 
is the velocity and £o a constant energy shift. Therefore, 
spin f is right moving (v k f = d k e k ^/h — vq > 0) and spin 
J, is left moving. The states ip k ^ and tp-ki form a Kramers 
pair, since Qipkf = +ifi-ki and <dtp-ki = -<Pkf, where 
G = —ia y K is the time-reversal operator consisting of a 
Pauli matrix a y and a complex conjugation operator K . 
It is possible to find a specific form of the HESs within the 
BHZ model, see Refs.l20l.l2ll and Appendix|B] Within this 
model, each of the two HESs in the Kramers pair consists 
of a mixture of two orbital states both with either positive 
or negative total angular momentum projection. Thus, 
it is possible to model the HESs as spin-1/2, which is 
sufficient for the present purpose. 

The coupling of the HESs to the spin bath is mod- 
eled as being point-like both along and transverse to the 
HESii^, i.e. V(x,y) = $ £. 6(x- X^y-Y^s-V , 
where J is the coupling (energy) constant, a is the area 
covered by a single spin S J = (S^, S^, S 3 Z ) at the fixed 
position (Xj,Yj) and s = (s x , s y , s z ) is the HES spin. 
Using the eigenstates in Eq.®, the interaction with the 
spin bath becomes 



k,k' j 



,i(k'-k)Xj 



i p qj J „ i p J „ 

T J k^,k'i°-k\ Lj k'i ' J kl,k'1- J + C kl c k't 



^itt,^'T^ C fet C fe't ^ki,k'l^z c ki C k'i 



where S± = S%. ± iS y are the raising and lowering oper- 
ators of the j th localized spin and ct (c fccr ) is the cre- 
ation (annihilation) operator of the HES (fka- The po- 
sition and wave-vector-dependent coupling is J 3 ka k , a , = 

^f/ S fe(^)/ f 'fe'(^)> where s = +1(-1) for a =t (|). 
The interaction consist of terms representing two kinds 
of scattering from k' to k of the j th localized spin: (i) 
Scattering within a HES leaving the spin bath unchanged 
(the two first terms) and (ii) scattering between the HESs 
by flipping a localized spin (the two last terms). It is the 
second kind of terms, which induce elastic scattering be- 
tween the HESs from k to —k. 

The factor e^ fe ~ k ) x i in Eq.© stems from the point- 
like nature of the interaction. If the spin-independent 
part of the potential was taken to be more extended in 
space, then the factor e^ fe ~ k ) x 3 would be replaced by 
the Fourier transform of the potential in the x direction 
at k' — k. Likewise a broadening in the y direction would 
introduce an integral over y in the coupling matrix ele- 
ments J kak i a i- However, these complications are not of 
importance for the basic physics of the elastic backscat- 
tcring between the HESs from k to —k discussed here, 
but could change their effectiveness. Furthermore, here 
an isotropic interaction between the HES spins s and the 
localized spins S J is used for simplicity. However, our 
results arc not affected, if uniaxial anisotrop y 12 i 13 i 25 is 
allowed [such that the coupling constants J would be 
different in the z direction and in the (x,y) plane]. 

In passing, we note that a small energy gap can open 
up in the HES spectrum due to the spin bath. A way 
to realize this, is by averaging over the positions of the 
localized spins in V, such that translational invariance is 
restored. Treating the sum of all the localized spins as a 
classical field, the Hamiltonian can be diagonalized and 
it becomes evident that the gap is basically proportional 
to the in-plane field. 



IV. INTER EDGE STATE SCATTERING RATES 

For weak coupling to the spin bath, the elastic inter- 
HES scattering rates T a4 _ a i can be found by the Fermi 
golden rule to be 

IVt = 2nh-± ]T \4 Uk ^\ 2 n k ^{l - n ki )5(e ki - e k ^), 



N- 



kk'j 



IV4 = 27rft-jy£ \ J kk'i\ 2 n k 'i(l - n kt )S(e k ^ - e k ^), 



kk' j 



(8) 



j J 



(7) 



where n ka is the electronic distribution function of the 
HES ct. A detailed derivation is given in Appendix 
The physical intuition behind these rates is clear: A scat- 
tering from, say, the HES f to 1 requires (i) an occupied 
state to scatter from, oc n k if, and (ii) an empty state to 
scatter into, oc 1 — n k \.. The scattering process conserves 
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energy, i.e. £14 = Ek'f- Moreover, the inter-HES spin- 
flip scattering t - M- demands the presence of a spin down 
in the bath to mediate the spin hip, which leads to the 
factor N±/N s . 

We observe that inserting these rates into the current 
Eq. © , the result is consistent with the Boltzmann equa- 
tion approach 2 ^— for general distributions n kcr . If the 
time it takes to move through the HES L/vq (the traver- 
sal time), is much shorter than the spin- flip time (inverse 
spin- flip rate) , then the distributions of the HESs n k( j are 
approximately equal to the lead distributions at which 
they originate, i.e. 



Magnetization versus bias voltage, r r = 



and 



riki * ft(e H ), (9) 



where the leads are Fermi distributed /°(e) = {1 + 
cxp[(e — fx a ")/k B T]} for a = L,R and T is the tem- 
perature. Therefore, using N^/N s = |(1 + M) and 
NJN S = |(1 - M), we end up with 

r T ^ = (l + M)r^, IV- t = (l-M)I*_ t , (10) 

where the rates at zero magnetization r^. < _ ., are 

J dk^ I Jl kiM \ 2 fl{e k ,)[l - /°(e fet )], 

(11a) 

(lib) 



hLv 



Here 



v = L/(2irhvo) 



the 



I P I 2 



density 
2 



of 



states 
4 



and 

[Ja/{2hL)Y[.f +k {Y ] )] i . Due 
to the energy window of the Fermi functions — /„], 
it is now explicitly clear that for M = scattering from 
the HES f to X dominates for [i L > [i R and vice versa as 
illustrated in Figs.QJb) and[TJc). 



V. THE MAGNETIZATION AND CURRENT 
CHANGE 

Now the steady state magnetization M s t is easily 
found by inserting the rates Eq. (|10|) into Eq. Q and solv- 
ing d t M = 0, i.e. 



1 1<- r 



1 t-f-4 



2 1 r 



(12) 



Furthermore, Eq.(|2]) gives that the magnetization builds 
up from being initially zero as M(t) = M s t{^ — e - */ r ™), 

1 is the character- 



where T„ 



7v s [2(r^ t +r^)+r r ] 

istic time scale for the magnetization process. Note that 
while magnetization builds up (i.e. in the non-stationary 
regime), t < r m , electrons are backscattered, causing a 
finite transient current change even for T r = — in con- 
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FIG. 2. (color online) The current-induced magnetization 
of the localized spins versus bias voltage over temperature, 
(fiR — fiL)/k B T, without additional spin-flip mechanisms in 
the bath, i.e. T r = 0. The simple approximation (|14|) for the 
magnetization (full black line) is compared to a numerical 
calculation using the transverse wavefunction from the BHZ 
model with N s = 10 3 randomly chosen spin positions Yj for 
£f = (fJ-L + I^r)/2 = and a temperature of 0.1K (purple 
stars), IK (red dots) and 10K (blue squares). 



trast to the stationary regime. 

Neglecting the weak additional spin-flip mechanisms in 
the bath, T r = 0, the magnetization is readily obtained 
from the rates T° a ^ a , in Eq. pTj) as 



M st = 



£jdfc|/ fc (y 3 )l 4 [/£( £fct ) 



Ej*l/fe(^)| 4 F( £fef ) 



(13) 



where F(e) = f° L (e)[l - /°(e)] + f° R (e)[l - was 
introduced. In the limit of bias voltage fiR — [1l and 
temperature k B T much smaller than the energy variation 
of the transverse eigenstate, the function |/fc(y))| 4 can 
be taken outside the integrals in Eq. ([T3")) , such that it 
simplifies to 



M s 



tanh 



{ V-L - m 
\ 2fc R T 



(14) 



Figure [5] shows that this is a very good approximation 
for a HgTc TI. Interestingly, Eo. (fl"4f resembles the well- 
known expression of a thermal equilibrium ensemble of 
spins in an external magnetic field^S, if the bias is ex- 
changed by the Zeeman energy. In contrast, Eci.([T4"]) de- 
scribes the current-induced magnetization at zero exter- 
nal magnetic field, i.e. a non-equilibrium steady-state sit- 
uation. From Eq. (TTJ| it follows that for bias much larger 
than temperature, maximal magnetization M s t — ±1 is 
found. Furthermore, in linear response — Hr\ *C fc B T, 
we have M st ^ (jjll - h t R )/(2k B T). 

To test the validity of the approximation (fH]l . it is 
compared to a numerical calculation of the magnetization 
Eq. (fT3)) in Fig. [2] To this end, we use the transverse state 
fk(y) from the BHZ model and randomly chosen spin 
positions Yj (see Appendix IBl for details on fk{y))- For 
the parameters for a 70A wide HgTe quantum well^, the 
magnetization expression (|14[) is found to be an excellent 
approximation for < T < 50K and N s from ten and up, 
see Fig. [5J Numerically, the magnetization is also found 
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to be independent of the Fermi level £p = + M_r)/2. more vanish for I\ = as expected. 

Next, we turn to the case of including a weak ad- 
ditional spin-flip mechanism, T r ^= 0. For increasing 
l/^fl — HL\/k B T beyond one, Eq. fTTj) clearly shows that 
one of the rates r° < _ CT , will increase while the other go 
to zero. Thus, for T r ^ 0, it is still possible to achieve 

M st a i (-1) for r°^ t > r r (r°^ 4 > r r ), which is 

a stronger requirement than — 3> k B T as in the 
r r = case Ea. fF!]) . Therefore, the maximal possible 
current change Eq.© is SI = ±eT r /2 for M s t — > ±1- 



The sample specific information about the positions 
(Xj,Yj) of the localized spins is not of importance. 
Therefore, we introduce the position average of a quan- 
tity A as A = ^n s J dXidYi ■ ■ ■ dXi\f B dYpf s A in ana- 
logue with impurity averaging 3 ^. In other words, it is 
assumed equally likely to find a localized spin j every- 
where in the area covered by the HESs. This enables 
us to give simpler expressions for magnetization and cur- 
rent change for r r 7^ 0. The position averaged inter-HES 
scattering rate for A4 = is 



'/ 



(15) 



where ns(e) = [e £ / k " T ~ 1] _1 is the Bose function and r\ = 
{■k/A)(uJ) 2 N 8 /N 2 is a dimensionless constant related to 
the strength of the interaction. Here N = LW y /a is the 
number of atoms covered in the plane by the HESs, since 
both an atom and a localized spin is taken to cover an 
area of size a. To obtain this, we used J °° dYj [/& (Yj )] 4 ~ 
1/Wy (see Appendix [B|l . For simplicity, the possibility of 
a weak k dependence of the width W y is neglected, which 
is justified within the BHZ model in Appendix [B] The 
opposite rate r5^_ , is found by interchanging [ij, and /j,r 
in F?^^ Eq. (|15|) . Therefore, using the position averaged 
rates, the magnetization Ea.(|12p becomes 



Mst = 



~ (J>R 



(M L -/x R )coth'^^ 



2fc„T 



(16) 



2)/ 



Thus, maximal magnetization is accessible for \hl — 
{Xr\ 3> 5^r r and — Mh| 3> k B T. The current change 
is readily found from Eq.([5]) to be 



SI = - — V 



TrhT r 



h e y coth (^ + Ar r ' 



(17) 



where the bias voltage V = (fiR — Mi)/ e was introduced 
such that l(°> = n-V, Therefore, in linear response the 
correction to the ballistic conductance — e 2 /h is 



SG 



irhTr 



h 2k B T + |r r : 



(18) 



i.e. both current and conductance are found to decrease 
compared to the case without a spin bath and furthcr- 



VI. DISCUSSION AND SUMMARY 

In this paper, we have discussed the current-induced 
magnetization at a single boundary of a 2D TI. For a two 
terminal device with two well-separated boundaries, the 
spin structure of the HESs is reversed on opposite bound- 
aries (i.e. if spin-up is right moving on the lower edge of 
the sample then it is left moving on the upper edge and 
vice versa). This means that the current- induced magne- 
tization has opposite signs - but equal magnitudes - on 
the two opposite boundaries. The bulk region remains 
non-magnetic and therefore, the sample as a whole is not 
magnetized for a two terminal symmetric setup (i.e. op- 
posite edges have the same length and spin concentra- 
tion). However, by adjusting the geometry, making the 
environmental spin concentration inhomogeneous or by 
using a multi-terminal sample an overall non-zero mag- 
netization can indeed be engineered. A possible appli- 
cation of the current-induced magnetization described in 
this paper, could be an all electrically accessible memory 
device. 

Two recent work o 32 i 33 on 3D TIs doped with magnetic 
impurities discuss how the surface of the TI can magne- 
tize, while the bulk remains non-magnetic (for a certain 
range of temperatures). In these studies, the magnetic 
ordering is caused by lowering the temperature leading to 
a phase transition in an equilibrium setup - in contrast 
to the present work describing current-induced magne- 
tization. Nevertheless, both effects build on the special 
nature of the edge and surface states in TIs, where the 
momentum and spin are locked together. 

Next, we estimate the magnetization time r m for mag- 
netic impuritie o 16 ! 17 to be on the order of tens of fj,s for 
T = 2K, (Xr — Hl ~ k B T, la — lcVA 2 and neglecting F r . 
Here the velocity vo ~ 4 x 10 5 m/s and width W y ~ 40nm 
correspond to a 70A wide HgTe quantum well within the 
BHZ modeli. Note that r m is independent of the HES 
length L and concentration of spin N s /N. If the local- 
ized spins are the nuclear spins of the 2D TI, then r m 
is much longer since J is much smaller. Therefore, the 
physics discussed in this paper is more relevant for mag- 
netic impurities embedded in a 2D TI. In both cases, 
the magnetization time is much longer than the traversal 
time: L/vq ~ lOps for L ~ 1/irn. Furthermore, these es- 
timates indicate that the traversal time is much shorter 
than the spin-flip time ~ r m /N s for low concentrations 
N s /N ~ 1CP 2 , such that inserting the lead distributions 
Eq.© into the rates Eq.© was indeed a justified ap- 
proximation. 

In summary, we have shown how a current through the 
HESs of a 2D TI can induce a magnetization of localized 
spins embedded in the TI. We have explained this by 
a simple physical picture of spin-flip scattering between 
the HESs. Only the region covered by the HESs mag- 
netize. We have demonstrated that if the spin bath is 
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only affected by the inter-HES scattering, then the sys- 
tem remains ballistic. In contrast, if an additional weak 
spin-flip mechanism is present in the bath, then a finite 
current decrease is found due to a competition between 
this mechanism and the spin-flip scattering between the 
HES. 
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Appendix A: Derivation of the inter helical edge 
state scattering rates 

In this Appendix, we give a detailed derivation of the 
clastic scattering rate Eq.© r CT ^ CT / of transferring an 
electron from the HES with spin a' to the HES with spin 
(7. In the limit of weak coupling to the spin bath, it can 
be found by the Fermi golden rule (see e.g. Ref. I31I). 

r — ' = Y \(f*\n<")\ 2 WiJ(E fa - E ial ). (Al) 

Here the initial state \i a ') (and final state \ f a )) consist of 
both an electronic part \i e ,a') and a part for the spin bath 
\is), i.e. \i a i) — \i e ,a') ® \is)- The Wi a , is the occupation 
factor of initial states^, and Ef a ~Ei a , is the total energy 
difference between the final and initial state. 

As a specific example, we calculate the rate r 4, <_-)-. To 
this end, we note that the final state is given in terms 
of the initial state as \f±) = c\> ±c kl ^.\i e rf) ® Ji^+\^s), i- c - 
an electron is transferred from k\ f to k[ 4- at the cost of 
a spin flip of the j' th localized spin. Therefore the sum 
over final states in Eq. (|Aip becomes a sum over ki, k[ 
and j'. The occupation factor is written as a product of 
the electronic and spin-bath part, Wi t = Wi c it Wi s , such 
that 

rVt = f E E \4i,k^\(is\S^Si\i 3 )\ 2 W is 

klk[ jj' 

fcfc' »e,t*S 

X K**e,tl4 1 t C *i; C U C fc'tl««.t)| aW *..t*( S /l -Bit)' ( A2 ) 

by inserting the interaction between the localized spins 
and the HESs Eq. (0 into the Fermi golden rule Eq. (|A1|) . 



The electronic part of the sum gives 

E 1 ( i e,ti4it c *i4. c U c *'ti i «.t> i 2w vt 

*<s,t 

= S ku k'S k ^ k n k ^(l - nki), (A3) 

where n ka is the electronic distribution function of the 
HES with spin a. The initial state of the spin bath is 
\is) = \mi, . . . , mjv s ), where irij = ±1/2 is the spin state 
of the j localized spin. Therefore writing the occupa- 
tion factor Wi s as a product over all the localized spins, 
W is = ■ ■ ■ W mNs , we obtain 

^Y.\^\ S ^ S i\^)\ 2 W ls =5^W mj ^ (A4) 

is 

where W mj= i is the probability that the j th spin is in the 
state rrij = —1/2. Here we used that the probabilities 
sum to one, i.e. W mi —\, + W mi= f = 1 for all i. We take 
the probability W mj= i to be W mj =i = N±/N s , where N a 
is the number of spins in the bath covered by the HESs 
with spin a and iV s = A|_ + JVf. Therefore, we end up 
with the rate 

kk'j 

The opposite scattering rate r^j. is found similarly and 
these are the results used in the main text in Eq.©. 



Appendix B: On the transverse eigenstate in the 
BHZ model 

A minimal model for describing the physics of a HgTe 
quantum well was proposed by Bernevig, Hughes and 
Zhang£. It is derived using k • p theory, see e.g. Ref. [l], 
ID, and H3 for further details. The BHZ model gives the 
following 4x4 Hamiltonian in 2D (k x , fc y )-space 

(Bl) 

in the basis {\E+),\H+),\E-),\H-)}. Here \E±) and 
\H±) arc Kramer pairs of electron-like and hole-like 
states, respectively. The states \E+) and \H+) (\E—) 
and \H— )) have positive (negative) angular momentum 
projection along the z direction, which is perpendicular 
to the 2D quantum well. The 2x2 block in H(k x , k y ) is 
given by 

^Ur-tit-^')' (B2) 

where e k = -D(k% + ky), M k = M - B{kl + k%) and A, 
D, Mo, B are parameters. From this model, Bernevig 
et a/.— predicted that HESs exist for HgTe quantum 
wells thicker than a certain critical value of 63A. Here 
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FIG. 3. (color online) (a). The energy dispersions eta 
(red/blue full lines) for the HESs in the BHZ model for a 
70A thick HgTe quantum well Eq.(B3)l. The HES exist in 
the bulk band gap E g — 2|Mo| and the bulk bands are 
shown in gray. (b). The transverse wave function Eq. (|B4a[) 
in the BHZ model for k = —O.OOlA -1 (brown dashed line) 
and k = — 0.005A -1 (purple full line), which have penetra- 



tion depths of A 2 



378 A and A" 1 



177A, respectively. 



The dotted (brown/purple) horizontal lines in (a) show the 
energy corresponding to these two transversal wave functions. 



we use the parameters for a quantum well thickness of 
70A (i.e. well within the TI regime )i: 



A = 3.65eVA, B = -68.6eVA 2 , 
D = -51.2eVA 2 and M n = -O.OleV. 



(B3) 



A priori the Hamiltonian in Ea. ljBll) has periodic 
boundary conditions and therefore does not contain any 
edges. However, by introducing boundaries in the model 
Hamiltonian (|B1|) . it is possible to find edge states in the 
TI regime. One way to do this is by introducing hard 
wall boundary conditions and using the Peierls substi- 
tution k a = —id a (for a = x, y) 20 ' 21 . The upper block 
h(k) leads to one HES, while the lower block h*(— k) 
leads to its time-reversed Kramer partner. Therefore, 
the HESs are a mixture of the orbital states with either 
positive {\E+} and or negative (\E—) and \H—)) 



total angular momentum projection. It is therefore that 
one can model the states as spin-1/2 in many situations, 
including the present one. In the limit of the bound- 
aries being very well separated, the dispersions become 
exactly linear—, £fc-Hn = £0 + ( — )^ok, where the en- 
ergy shift and velocity are given by Eq = —MqD/B and 
v = -V-B 2 - D 2 \A\/(HB), respectively, see Fig. E^a). 
For a hard wall boundary at y = (and the TI at y > 0), 
the (real and normalized) transverse wave function is 
found to be2&2i 



) , (B4a) 



(Ai - A 2 ) 2 

where the k dependence is in the parameters Ai and A 2 : 

(B4b) 
(B4c) 




Here we introduced 
4 B 



(B 2 - D 2 ) 



D\A\VB 2 - D 2 
B 



k + (B 2 - D 2 )k 2 



(B4d) 



Note that the transverse wave function vanishes at the 
boundary, fk(0) = 0, as required. It is evident from the 
form of the transverse state that the penetration depth or 
width is given by W y = A^ 1 (since A2 < Ai). Figure Mb) 
shows this transverse state for the two different values of 
the wave vector k corresponding to an energy difference 
as large as about half an energy band gap (~ lOmeV). 

In the main text, we use the transverse wave func- 
tion with the parameters in Eq. ([B3[) as an example. We 
use it for numerical evaluation of the magnetization in 
Fig. [2] of the main text. Furthermore, in the averag- 
ing over the positions of the localized spins, we use that 
J °° dYj[fk(Yj)] 4 ~ 1/Wy. This can be justified by using 



the BHZ transverse wave function Ea. (|B4aj ) such that 

3AiA2(Ai + A2) 



dY j \fk<Y,W = 







(3Ai + A 2 )(Ai +3A 2 ) : 



(B5) 



which in the limit A 2 

14 ~ a 



<C Ai gives the used result 
dYjlfktYj)]* ~ A 2 = 1/Wy. Moreover, the weak k 
dependence of W y is neglected in the main text, which is 
seen to be reasonable both from the explicit form (|B4|) 
and the examples given in Fig. EHb). 
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